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Abstract. - The influence of dissipation on the fluctuation statistics of the total energy is inves- 
tigated through both a phenomenological and a stochastic model for dissipative energy-transfer 
through a cascade of states. In equilibrium the states obey equipartition and the total energy obeys 
the central limit theorem, giving Gaussian fluctuation. In the presence of dissipation the fluctua- 
tions can be driven non-Gaussian if there is macroscopic energy transfer from large to small scales. 
We are thus able to equate the non-Gaussian order parameter fluctuations in model equilibrium 
systems at criticality with energy fluctuations in these dissipative systems. Energy fluctuations in 
the phenomenological model map directly onto the l//"-noise problem and numerical simulations 
of the stochastic model yield results in qualitative agreement with these predictions. 
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Introduction. — Dissipative many body systems such 
as turbulent flows, granular matter or evolving dry foams 
Qan be parameterized by spatially averaged energy func- 
'ions for the integrated kinetic energy, energy dissipation, 
qr injection in the case of driven systems. Despite the 
many body nature of these systems, the probability den- 
sity for such global quantities is often found to be non- 
Gaussian (for example see [IHI3]), with finite skewness 
4nd a greater probability for large deviations than pre- 
dicted by the central limit theorem. Although the skew- 
iiess varies from one system to another, a generic pattern 
does seem to emerge |14| in which the energy distribu- 
tion can often be approximated by a generalized Gumbcl 
function [15H18j -a distribution usually appearing in the 
context of extreme value statistics [l9l[20] . Such functions 
also parameterize critical order parameter fluctuations in 
a class of Gaussian models sitting at their lower critical 
dimension. In particular they give the form of order pa- 
rameter fluctuations along the line of critical points of the 
2D-XY model to an excellent approximation [T4HI6] . as 
well as the form of one-dimensional interface fluctuations 
showing 1// noise spectra [3T1[25]. In these equilibrium 
systems, while the probability for order parameter fluctu- 
ations is non-Gaussian, energy fluctuations perfectly sat- 
isfy the central limit theorem, as the energy statistics for 
individual elements is given by equipartition. Moving the 
no n- Gaussian fluctuations from order parameter to energy 
suggests a breakdown of equipartition and the absence of 



a unique effective temperature covering all scales of the 
system. This is the subject of this paper. In the next 
section, using a simple model of a dissipative system, with 
energy transfer between scales and scale dependent dis- 
sipation, we show how dissipation induces a non-uniform 
energy profile among scales and drives the system towards 
a steady state with scale dependent temperature. In the 
third section, we test these ideas against a simple stochas- 
tic cascade model which we simulate numerically. Finally, 
in the discussion, we argue that for a dissipative systems 
with macroscopic energy transfer, the concept of equipar- 
tition of energy at equilibrium should be replaced by an 
equipartition of energy transfer in the steady state. 

A generic phenomenological model. In the fol- 
lowing we construct a simple model which develops the 
above phenomenology of non-Gaussian energy fluctuations 
and spatial correlations when driven out of equilibrium. 
The model consists of a large number of Fourier modes (j)q^ 
with q a wave vector, such that the energy associated with 
mode q is Eg = ^|(/)gp, considered as a random variable. 
The modes with different q values are considered to be 
statistically independent. In the spirit of a local (in scale) 
equilibrium approximation, we consider that all distribu- 
tions Pq{sq) have an equilibrium form, with a q-dependeiit 
temperature Tg, 
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a scenario which has been shown to be exact in a cas- 
cade model similar in spirit to the one considered here 
|23) . At equilibrium, one has Tq = Tq for all q, so that 
the total energy E = Eq is simply a sum of indepen- 
dent and identically distributed variables. From the cen- 
tral limit theorem, it follows that the distribution of the 
total energy is Gaussian in the limit of a large number of 
modes. Dissipative fluid systems, such as turbulent flows 
and granular gases can be modeled at this level by consid- 
ering an energy cascade process from large to small length 
scale, with dissipation occurring principally at small scales 
so that the cascade terminates at the dissipative scale qq. 
Note that phenomenologically, a one-dimensional cascade 
is often considered. In this case however, one should take 
into account the fact that the modes are space filling in 
d dimensions, so that each level of the cascade consists of 
a shell of order g'^"^ modes 0q. We shall come back to 
this point in the discussion. For the present, we do not 
take this density of state issue into account, and consider a 
purely one-dimensional cascade. For simplicity, we assume 
that energy transfer is diffusive, driven by a temperature 
gradient between scales, so that the diffusive energy flux 
is given by 

Jdif(g) = -Diq)^ (2) 

where D{q) is a scale-dependent diffusion coefficient. Note 
that we have implicitly replaced the discrete set of tem- 
peratures Tq by a function T{q) of the continuous variable 
q defined on the interval 1 < q < qo- In order to model a 
scale invariant system, it is natural to choose a power-law 
form, 

Diq)^Doq'+''. (3) 

For a ~ 0, the diffusive flux reduces to Jdii{q) = 
—dT/dhiq, so that the dynamics is exactly the same at 
all scales. The parameter a thus quantifies the dynamical 
bias between large and small scales. Injection takes place 
at g = 1, and we further assume that dissipation is local- 
ized at (7 = go 3> 1, so that energy is locally conserved in 
the interval 1 < q < qo. The conservation equation reads 



dT 
'dt 



d 



(4) 



leading to the following diffusive evolution equation for 
the temperature T{q): 



dT a ( dT 



(5) 



This equation is supplemented by the boundary conditions 
at q = 1 and at g = (jfo- At g = 1 the injection mechanism 
is modeled by assuming that it fixes the average value of 
the large scale energy: T{q = 1) = To. Ai q = qo, the flux 
of dissipated energy 



Jdis = i^T{qo), 



(6) 



is equal to the diffusive flux Jd\i{qo)- The dissipation co- 
efficient V is assumed to be much larger than the diffusion 



coefficient D{qQ), consistently with the assumption that 
dissipation is localized at the scale go- 

We now focus on the stationary state of the model; 
dT/dt = 0, so that Jdif(g) = Jo, independently of scale. 
This condition enforces a temperature gradient between 
scales in the out-of-cquilibrium steady state, thus ensuring 
that the system cannot support a single, effective temper- 
ature when driven away from equilibrium. It follows that 



dT 



Jo 



Doq^ 



(7) 



Taking into account the boundary condition T{q = 1) = 
To, the temperature profile is given after integration by 



T{q) - To 



Jo 
aDr 



Jo 
aDo 



(8) 



Using the equality of diffusive and dissipative fluxes at 
q = qo, we can determine Jo, yielding 



— To 



aDo 



(9) 



The resulting value of the flux thus depends on a. If a > 
0, given that qo ^ 1 and Dq/v 1, a finite flux Jo is 
obtained, 

Jo « aDoTo, (10) 



while for a < 0, the flux becomes very small 

\a\DoTo 



Jo 



(11) 



This transition between two regimes as a function of the 
parameter a is similar to that reported in f23| . From 
Eq. one finds that the temperature profile T{q) be- 
haves for a > as a power law of q. 



T{q) 



To 
9"' 



In contrast, for a < 0, the temperature profile is 



T{q) « To 



1 - 



(12) 



(13) 



and is thus essentially constant, and equal to Tq, unless 
q is close to go (on a logarithmic scale). Note that these 
asymptotic results rely on gj,"' being large. For small a, 
one thus expects strong finite size effects, as very large 
values of go may be needed to fulfill gj,"' » 1. The tem- 
perature gradient is therefore seen to be intimately related 
to finite energy transfer down the scales and hence to en- 
ergy dissipation at small scale. For negative values of a, 
although the transfer becomes small for large go, it remains 
sufficient to establish a local equilibrium on a sufficiently 
long time scale. 

We now turn to the fluctuations of the total energy 
E = J2q^q^ where is distributed according to Eq. ((T|), 
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Fig. 1: Illustration of the 1//" -noise distributions for several 
values of a. Inset: same distributions on a semi-logarithmic 
scale. 



with the temperature Tq satisfying either Eq. or 
Eq. ([T^ . It is clear that the energy E is a sum of indepen- 
dent random variables with non-identical distributions. If 
a < 0, the distributions of Eq are almost the same for 
all q, so that the central limit theorem holds in the limit 
qo — > cxD. In the opposite case a > 0, the distribution of 



£g reads 



Pq{eq) = nq''e-^'^ (14) 

with K = I/Tq, and one recovers the so-called l//"-noise 
problem investigated in [35]. To determine the shape of 
the energy distribution, wc introduce the rcscalcd energy 
X defined as 



(15) 



where (E) and as are the empirical mean values and stan- 
dard deviations of the energy E in each case. 

One might expect non-Gaussian statistics to appear for 
all a > 0. However, it turns out that the distribution P{x) 
remains Gaussian in the large N limit as long as a < ^ 
|22| . For a > i, the limit distribution is non-Gaussian, 
and becomes more and more asymmetric with increasing 
a, crossing over from Gaussian for a ^ to the expo- 
nential distribution for a — >■ oo. For a = 1, one obtains 
the Gumbel distribution [5T1[53], originally known in the 
context of extreme value statistics [Tni[lD|. Remarkably, 
it seems that experimental and numerical systems often 
present results that can be described to a good approxi- 
mation by a generalized Gumbel distribution close to this 
form [IHSI- The no n- Gaussian distributions are illustrated 
in Fig. [1] From the results of [12] , one easily obtains the 
skewness and kurtosis of these distributions as 



(x4)-3 = 



2C(3a) 
C(2a)3/2 

6C(4a) 
C(2^ 



(16) 
(17) 



where ^ is the Riemann Zeta function, defined as ^(a) = 
Note that by definition (x) = and (x^) = 



1. Both the flatness and the kurtosis go to zero when 
a ^ i"*", consistently with the fact that the 1 //"-noise 
distributions converge to the Gaussian law in this limit. 

To sum up, the main outcome of the generic phenomeno- 
logical model considered in this section is a transition be- 
tween Gaussian and non-Gaussian energy distributions, 
occuring when the parameter a characterizing the trans- 
fer properties between scales is varied. For a finite number 
of modes, the transition should be replaced by a crossover 
between close-to-Gaussian and strongly non-Gaussian dis- 
tributions. It is plausible that the precise crossover value 
of a may depend on the details of the model, but one 
may expect the scenario to be qualitatively valid in more 
generic situations. In the next section, we test this sce- 
nario in an explicit stochastic model. 

Stochastic dissipative model. — As a test of the 
above scenario, wc consider a dissipative energy transfer 
model explicitly defined by stochastic rules. This model 
is similar in spirit to the ones considered in [T7l [231125] . 
though some significant differences are present. The ener- 
gies are positive real variables Sq attached to the wavenum- 
bers q ~ l,...,(7max, forming a one-dimensional lattice 
with an energy cascade from low to high wavenumber. As 
in the phenomenological model, the energy of mode q is 
Eq = where (j)q is a Fourier amplitude. The dy- 

namics is defined by three different mechanisms: energy 
injection at large scale, internal energy transfer between 
neighboring scales and local energy dissipation. To pre- 
serve scale invariance, the presence of energy transfer be- 
tween wavenumbers q and q' is determined according to 
the ratio q'/q; namely it occurs only if 1/A < q'/q < A, 
where A > 1 is a fixed parameter. The different trans- 
fer rates are defined as follows: an amount of energy ^ 
is either injected at q = 1 with rate (probability per unit 
time) I{fJ-), transferred from a site q to site q' with a rate 
Tq,qi {fi\eq), OT dlssipatcd from site g at a rate Aq{fi\sq). 
Note that the transfer and dissipation rates depend on 
the energy Sq present at wavenumber q. Wc choose the 
following forms for the rates 



Tqy{n\eq) 
Aq{^l\eq) 



Doiqq') 



'1 ' 



(18) 
(19) 
(20) 



Choosing this form for /(^) means that energy is injected 
at wavenumber g = 1 by connecting this element to a 
thermostat at temperature Tq = 1, which fixes the energy 
scale. Energy is also transferred back to the reservoir with 
a rate I{^\sq) = Sq. The value of the parameter a is again 
the key to the way in which energy is transfered between 
scales. The constant Dq in Eq. (fT9]) is chosen such that 
Doiqq')"^^ remains bounded when (/max oo. We take 
Z^o = 1 for a < and Dq = g~ax for a > 0. Finally, 
the constant vq controls the dissipation rate and can be 
interpreted as a viscosity. In the following, we shall be 
interested in the limit where vq takes small, but non-zero 
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values. Note that the dissipation rate is the only dynam- 
ical rule that breaks detailed balance. In the absence of 
dissipation (i^o = 0) the system converges to the equilib- 
rium Gibbs distribution with temperature Tq = Tq at all 
scales. We shall sec however that the limit f o —5" is sin- 
gular, as is well-known in the context of fluid turbulence 

The main differences between the present cascade model 
and that introduced in [17j are that in the latter the en- 
ergy transfer rate is fully biased and scale-independent 
and that the dissipation rate has an ad hoc form chosen 
to make the model analytically solvable. Here, in addition 
to the inclusion of back transfer from large to small scale, 
we choose a dissipation rate proportional to g^, which is 
natural in the context of viscous damping. 

With this form, as long as a < 2, the ratio of the dissi- 
pated to transferred energy increases as the wavenumber q 
increases, assuring a cut-off wavenumber for the cascade: 



90 



l/(a-2) j^a/(a-2) 



(21) 



In the context of simulations of a finite system with a 
number gmax of modes, the characteristics of the cascade 
are determined by both qo and (/max and their interplay, 
as we show below. 

We have simulated the model for gmax — 500, and mea- 
sured the steady-state distribution of the total energy for 
several values of a. The resulting distributions P{x) of the 
rescaled energy x are shown for i^q = 10~^ in the upper 
panel Fig. [21 on a linear scale in the main figure and log- 
linear scale in the inset. For a = 2 the shape is strongly 
non-Gaussian, with a clear skewness and long tail for fluc- 
tuations for positive values of x. As a decreases from 2 
to 1 there is a sharp decrease in the anisotropy of the dis- 
tribution in qualitative agreement with the data for the 
phenomenological model shown in Fig. [T] The form for 
a = 1 is close to that typically observed in experiments 
and simulations for dissipative systems (for example see 
[THH] ) . Further reduction in the value of a sees a slow evo- 
lution towards a Gaussian distribution although it is not 
possible to determine clearly from this data if a change in 
regime occurs for a around zero. 

Evidence for such a crossover is however given by the 
evolution of the energy spectrum (e^) with a, shown in 
the lower panel of Fig. [5] There is a clear quantitative 
difference between the spectra for positive and negative 
values of a. For a = 2, the spectrum is approximately 
power law over five orders of magnitude in energy (with 
an exponent a w —2.9), with no evidence of a crossover 
to a dissipative regime, which is consistent with the ratio 
of energy transfer to dissipation being scale independent. 
For a = 1 the inertial range of approximate power law 
scaling (here with an exponent a « —1.6) is cut off by a 
build up of energy at a scale around q = 2 x 10^, which 
corresponds to the dissipation scale, go for this value of 
vq- One might have expected a fall off in energy per mode 
as one goes to larger wavenumber. However, in this case 




Fig. 2: Top panel: Distribution of the rescaled energy x — {E — 
{E))/(te in the stochastic cascade model for different values 
of the a characterizing energy transfer through scales (vq = 
10"'^;. Inset: same data on a semi-logarithmic scale. The 
dotted line is the Gaussian distribution, shown for comparison. 
Bottom panel: energy spectrum (e,) for the same values of a, 
on a log-log scale. Straight lines indicate slopes of —2.9 (a = 2) 
and —1.6 (a = 1). 



go is close to g,nax with the result that energy is stored at 
the bottom of the cascade, rather than being dissipated 
at larger wavenumber. Note that when simulations are 
run for smaller values of vq, we find that the measured 
exponents a decrease in magnitude (data not shown) to- 
wards a, as for the phenomenological model presented in 
the previous section. More extensive finite size scaling is 
clearly necessary here in order to capture the details of 
the cascade and make a quantitative comparison between 
the phenomenological and the numerical model, but the 
general behavior of extensive energy transfer, inertial scal- 
ing regime and non-Gaussian energy fluctuations appears 
well established. For negative values of a the spectrum 
again mirrors that for the phenomenological model, being 
constant at small q and cut off to zero around the charac- 
teristic dissipation scale, go ^ 50 for a = — 1 and go ^ 15 
for a = -2. 

The deviation from the Gaussian distribution for a < 
can be traced to flnite viscosity effects. We have computed 
the distribution P{x) for decreasing values of the viscosity 
for a = —2, shown in Fig. [3] As in this case we 



estimate the crossover scale go 



-1/4 



, we should be able 



to push vq as low as 10 ^ before the maximum cut off 
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Fig. 3: Distribution P{x) of the rescaled energy x for decreasing 
values of the viscosity uq (a — —2). A slow convergence to the 
Gaussian distribution (dotted line) is observed. Inset: average 
energy (e,) as a function of q for the same values of uq (vo 
decreases from bottom to top). 



'/max = 500 plays a significant role. A slow convergence to 
the Gaussian distribution is observed, as I'o decreases. The 
convergence is indeed associated with the increase of the 
caracteristic decay scale of the energy spectrum {sq) when 
fo — > 0, go (see inset of Fig. [3]), which does indeed scale 
roughly as multiplied by a factor of order unity. This 

evolution is consistent with an increase in the number of 
degrees of freedom which are equilibrated at temperature 
To, allowing an approach to Gaussian behavior for P{x) 
as goes to zero. 

To provide further evidence for a crossover from Gaus- 
sian to non-Gaussian behavior as a function of a, we com- 
pute the skewness (x^) and the kurtosis (a;*) — 3, as a 
function of a for fixed = 10^^. Results are plotted in 
the Fig. m The skewness and kurtosis take small values 
for a < 0.5, and larger values for a > 0.5. The data 
can be compared with the analytic predictions from the 
phcnomcnological model. Although we are not able to 
make a quantitative statement, the simulation results are 
clearly consistent with the phenomenology within the ac- 
curacy achieved. For a < 0, the skewness and kurtosis 
can be made smaller by decreasing the viscosity vq, as 
shown in Fig. [31 In contrast, for positive values of a, 
numerical simulations are difficult for very low viscosities, 
as gmax must exceed go to avoid energy accumulation at 
large wavenumber. A more detailed study therefore re- 
quires finite size scaling both in viscosity and gmax which 
we leave for a future study. However the results presented 
in this section show that the scenario of an energy trans- 
fer driven crossover between near-Gaussian and strongly 
non-Gaussian energy fluctuations is also present in this 
stochastic cascade model. 

Discussion. — In the previous sections we have shown 
that the asymmetry of the energy distribution in dissipa- 
tive systems can be traced back to the efficiency of the 
transfer mechanism from large to small scales. This de- 
parture from the central limit theorem has, in the past 



Fig. 4: Skewness and kurtosis as a function of a for the sim- 
ulated model (symbols) and for the l/f"' -noise distributions 
(lines), with a viscosity vo = 10~^. 



been compared with non-Gaussian order parameter fluc- 
tuations in model systems of criticality at equilibrium 
[T4lll6l[2T] . The asymmetry can here be understood as re- 
sulting from the constraint imposed by the equipartition 
of energy. Consider the Fourier amplitude of the order 
parameter for wave vector q, ?/;q in the specific example of 
the low temperature phase of the 2D-XY model: the en- 
ergy associated with the mode reads eq = ^ g^li/'qP, while 
the contribution of wavevector q to the order parameter 

is TTlq = i IV'qP 



|16] . The equipartition of energy at equi- 
librium imposes that the average energy (eq) of mode q 

fcsT. The equipar- 
so that the scale of 



is independent of q, namely (cq) = 
tition thus imposes (wq) = fc^T/g^ 



contributions to the order parameter over the Brillouin 
zone diverges in the thermodynamic limit. Combining the 
1/g^ scaling for mq with the density of states, g{q) ^ q 
results in non-Gaussian fluctuations for the order parame- 
ter, while the fluctuations of the total energy are perfectly 
Gaussian, being the sum of independent and identically 
distributed random variables. Similar behavior occurs in 
one dimension if interactions are long ranged |27] , leading 
directly to a " 1/ /-noise" spectrum [5T] , which can in turn 
be generalized to the 1//" behaviour [55] discussed above. 
This is of course an extremely simplified vision of critical- 
ity |28| : critical systems are certainly not in general inte- 
grable and non-linear coupling between modes can drive 
a divergent specific heat and consequently non-Gaussian 
energy fluctuations. However, the specific heat exponent 
a is systematically smaller than that for the susceptibil- 
ity, 7, indicating a weaker divergence and showing that 
energy fluctuations are in general more constrained than 
those for the order parameter. 

Here we consider dissipative systems at the same level of 
approximation, finding clear violation of the central limit 
theorem and non-Gaussian statistics for the spatially aver- 
aged internal energy. In the context of an energy cascade, 
with dissipation localized at small length scales, equipar- 
tition is replaced, in steady state, by a new constraint; the 
average energy transfer between scales is uniform through- 
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out the cascade, Jq = Jq. For the models with diffusive 
dynamics that we have considered, this steady state con- 
dition ensures that global equipartitiou of energy must be 
violated, thus introducing a scale dependent temperature, 
Tq ~ and non-Gaussian energy fluctuations, if a is 
large enough. Indeed, as = one finds for both the 

phcnomcnological and stochastic models that Jq ^ q°'\(f)q\'^ 
in analogy with the energy of the equilibrium system. 
Hence we find quantitative equivalence between order pa- 
rameter fluctuations in the equilibrium models and those 
for internal energy in the case of dissipation. The sum 
over all instantaneous energy transfers would thus satisfy 
the central limit theorem and this rather non-standard 
global quantity would have Gaussian statistics. There is 
once again a barrage of simplifications separating our phe- 
nomenology from real systems, but our arguments do lay 
out a general framework for why energy fluctuations are 
more violent in dissipative systems than in equilibrium 
and in what context one can expect energy fluctuations to 
be approximated by l//"-like distributions. In particular 
examples, such as kinetic energy fluctuations of a granu- 
lar gas in a homogeneous cooling state [TU|, or of decay- 
ing Burgers turbulence [7] , the agreement between simula- 
tion results and the model distributions is close enough to 
suggest that the phenomenology presented here becomes 
quantitatively correct, at least in these specific cases. 

The comparison with equilibrium systems can be 
pushed further by considering the role of spatial dimen- 
sion, d in the presence of dissipation. At a heuristic level, 
a naive reasoning involves treating modes at each wave 
vector as statistically independent and consists in simply 
multiplying the 1//" spectrum by the density of states 
g[q) ^ q'^^^, yielding an effective exponent a' = a—{d—l) 
for the energy spectrum. If a is large enough, a transition 
from non-Gaussian to Gaussian fluctuations is observed 
when increasing d beyond d^ = a + (corresponding to 
a' = ^), which can be interpreted as an upper critical di- 
mension, in qualitative analogy with equilibrium critical 
phenomena. Taking into account correlations between the 
mode at each level of the cascade will however reduce the 
number of effective degree of freedom, probably placing 
the effective exponent in the range, a — {d — 1) < a' < a. 

In conclusion, our main results are two-fold. Firstly, 
we have shown that there generically exists a crossover 
from Gaussian to strongly non-Gaussian fiuctuations of 
the spatially averaged energy, when varying the transfer 
properties through scales (quantified here by the param- 
eter a) in dissipative systems. Secondly, we have argued 
that in the non-Gaussian regime, the asymmetry of the 
energy distribution is driven by a finite energy flux cross- 
ing the system from injection to dissipative scales. The 
constancy of this flux across scales plays a role similar to 
that of energy equipartitiou. In consequence non-Gaussian 
energy fluctuations arise in complete analogy to those for 
the order parameter at equilibrium. As for future work, it 
would be of interest to try to quantitatively validate this 
scenario in more realistic models (shell models. Burgers 



turbulence, granular gases...). 
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